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Precise calculations of core properties in heavy-atom systems which are described by the opera- 
tors heavily concentrated in atomic cores, like to hyperfine structure and P,T-parity nonconservation 
effects, usually require accounting for relativistic effects. Unfortunately, completely relativistic treat- 
ment of molecules containing heavy elements is very consuming already at the stages of calculation 
and transformation of two-electron integrals with a basis set of four-component spinors. 

In turn, the relativistic effective core potential (RECP) calculations of valence (spectroscopic, 
chemical etc.) properties of molecules are very popular because the RECP method allows one to 
treat quite satisfactory the correlation and relativistic effects for the valence electrons of a molecule 
and to reduce significantly the computational efforts. The valence molecular spinors are usually 
smoothed in atomic cores and, as a result, direct calculation of electronic densities near heavy nuclei 
^ | is impossible. 

O ' In the paper, the methods of nonvariational and variational one-center restoration of correct shapes 

of four-component spinors in atomic cores after a two-component RECP calculation of a molecule 

(-H , are discussed. Their efficiency is illustrated in correlation calculations of hyperfine structure and 

parity nonconservation effects in heavy-atom molecules YbF, BaF, TIF, and PbO. 

SHORT NAME: Studying core properties in molecules 
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INTRODUCTION 



Historical background of study P- and T-parity nonconservation effects in heavy-atom molecules. After discovery 
of the combined charge (C) and space (P) parity violation, or CP-violation, in if°-meson decay the search for the 
. electric dipole moments (EDMs) of elementary particles has become one of the most fundamental problems in physics 
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Permanent EDM is induced by the weak interaction that breaks both the space symmetry inversion 
and time-reversal invariance (T) ,.6| . Considerable experimental effort has been invested in measuring atomic EDMs 
induced by the EDMs of proton, neutron and electron and by P,T-odd interactions between them. The best available 
restriction for the electron EDM, d e , was obtained in the atomic Tl experiment 7], which established an upper limit 
of \d e \ < 1.6 x 10 -27 e-cm, where e is the charge of the electron. The benchmark upper limit on a nucleus EDM is 
obtained in atomic experiment on 199 Hg |djjg| < 2.1 x 10~ 28 e-cm, from which the best restriction on the proton 
EDM, \d p \ < 5.4 x 10~ 24 e-cm, was also recently obtained by Dmitriev & Sen'kov 0. 

Since 1967, when Sandars suggested to use polar heavy-atom molecules in the experimental search for the proton 
the molecules are considered as the most promising objects for such experiments. Sandars also noticed 
that the P- and P,T-parity nonconservation (PNC) effects are strongly enhanced in heavy-atom systems 
due to relativistic and other effects. For example, in paramagnetic atoms the enhancement factor for an electron 
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EDM, d a tom/d e , is roughly proportional to a 2 Z s a.D, where a « 1/137 is the fine structure constant, Z is the nuclear 
charge and an is the atomic polarisability. It can be of order 100 or greater for highly polarizable heavy atoms 
(Z > 50). Furthermore, the effective inner molecular electric field E mo \ acting on electrons in polar molecules can be 
a few orders of magnitude higher than the maximal field E cyLt accessible in a laboratory, E mo i/E ext ~ 10 5 . The first 
molecular EDM experiment (on search for the proton EDM and other nuclear P,T-odd effects) was performed on TIF 
by Sandars et al. [lJ (Oxford, UK). In 1991, in the last series of the 205 T1F experiments by Hinds et al. (Yale, 
USA), the restriction d p = (— 4± 6)xl0~ 23 e-cm was obtained (that was recalculated in 2002 by Petrov et al. |bj] as 
d p = (-1.7±2.8)xl0- 23 e-cm). 

In 1978 the experimental investigation of the electron EDM and other PNC effects was further stimulated by 
Labzowsky et al. [15L llr^ and Sushkov & Flambaum wn0 clarified the possibilities of additional enhancement 
of these effects in diatomic radicals like BiS and PbF due to the closeness of levels of opposite parity in fi-doublets 
having the 2 n 1 / 2 ground state. Then Sushkov et al. Q| and Flambaum & Khriplovich Q| have suggested to use re- 
doubling in diatomic radicals with the 2 Si/2 ground state for such experiments and the HgF, HgH and BaF molecules 
were first studied semiempirically by Kozlov 20] . At the same time, the first two-step ab initio calculation of PNC 
effects in PbF initiated by Labzowsky was finished in St. -Petersburg (Russia) [2]]]. A few years later, Hinds started 
experimental search for the electron EDM on the YbF molecule, on which the first result was obtained by his group in 
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2002 (Sussex, UK) 22], d e = (— 0.2 ± 3.2) x 10~ 26 e-cm. Though that restriction is worse than the best available now d e 
datum calculated from the Tl experiment (see above), nevertheless, it is limited by only counting statistics, as Hinds 
et al. pointed out in [22]] . 

New series of the electron EDM experiments on YbF by Hinds group (London, UK) are in progress and new 
generation of the electron EDM experiments using a vapor cell, on the metastable a(l) state of PbO, is prepared 
by group of DeMille (Yale, USA). The unique suitability of PbO for searching the elusive d e is demonstrated by the 
very high statistical sensitivity of the Yale experiment to the electron EDM. In prospect, it allows one to detect d e 
of order of 10 -31 e-cm 23], three-four orders of magnitude lower than the current limit quoted above. Some other 
candidates for the EDM experiments, HgH, HgF, and TeO*, are yet discussed and the experiment on PbF is planned 
in Oklahoma, USA. 

In order to interpret the data measured in such molecular experiments in terms of the electron EDM or other 
fundamental constants of P- and P,T-odd interactions, high-precision calculation of the electronic wave function near 
a heavy nucleus is required. Moreover, ab initio calculations of some molecular properties are usually required already 
prior to the stage of preparation of the experimental setup. 

Heavy-atom molecules, computational strategies. The most straightforward way for calculation of electronic struc- 
ture of molecules containing heavy atoms is solution of the eigenvalue problem using the Dirac-Coulomb or Dirac- 
Coulomb-Breit Hamiltonian |2 ll. l25l l2q] when some approximation for the four-component wave function is chosen. 

However, even applying the relativistic SCF approximation to heavy-atom molecules when solving the Dirac-Fock 
or Dirac-Fock-Breit equations followed by the transformation of the two-electron integrals is not always an easy 
task because much larger sets of primitive basis functions are required [27I for such all-electron four-component 
calculations as compared to the nonrelativistic case. Starting from the Pauli approximation and Foldy & Wouthuysen 
transformation, many different two-component approaches were developed in which only large components are treated 
explicitly (e.g., see |29j, yOj, |3l|, |32( and references). 

The most computational savings are yet achieved when the two-component relativistic effective core potential 
(RECP) approximation proposed by Lee et al. [3^] is used [3 13- There are the following reasons for it. The RECP 
approach allows one to exclude large number of chemically inactive electrons from molecular calculations and to treat 
explicitly the valence electrons only. Then, the oscillations of the valence spinors are usually smoothed in heavy-atom 
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cores. As a result, the number of primitive basis functions can be reduced dramatically. This approach is based on a 
well-developed nonrelativistic technique of calculations, however, an effective spin-orbit interaction and other scalar- 
relativistic effects arc taken into account usually by means of radially-local or separable potentials 

Many 

complications of the Dirac-Coulomb(-Brcit) molecular calculations |37j are avoided when employing RECPs. The 
radially-local RECP approaches for "shape-consistent" (or "norm-conserving") pseudoorbitals are the most widely 
applied in calculations of molecules with heavy elements though "energy-adjusted/consistent" pseudopotentials by 
Stuttgart team and Huzinaga-type "ab initio model potentials" 0| are also actively used. (In plane wave calculations 
of many-atomic systems and in molecular dynamics, the separable pseudopotentials |.3(1 l39l l4Cl l4ll |42( are more 
popular now because they provide linear scaling of computational effort on the basis set size.) The nonrelativistic 
shape-consistent effective core potential has been first proposed by Durand & Barthelat and then a modified 
scheme of the pseudoorbital constructing was suggested by Christiansen et al. and by Hamann et al. 0. 

However, inaccuracy of the conventional radially-local RECP approaches reaches 1000-3000 cm -1 for the transition 
and dissociation energies that may be insufficient in practice. In deep reorganization of electronic configuration 
structure of molecules containing, in particular, transition metals, lanthanides, actinides, and superheavy elements 
the conventional radially-local RECP approaches can not provide reliable results for a wide variety of properties (as 
it was shown both in many calculations and theoretically, see and references) though otherwise is sometime 

stated (e. g., s ee l48l l49l ) . Such problems can be efficiently overcomed by applying the generalized RECP (GRECP) 
approach \iA It?! IfiflflsiT , that involves both radially- local, separable and Huzinaga-type [ssL I53I 0] potentials as 
its components and as particular cases. In the GRECP concept, the inner core, outer core and valence electrons are 
first treated employing different approximations for each. 

Nevertheless, calculation of such properties as electronic densities near nuclei, hyperfine structure constants, P,T- 
parity nonconservation effects, chemical and isotopic shifts etc. with the help of the two-component pseudospinors 
smoothed in cores is impossible. For evaluation the matrix elements of the operators singular on nuclei, proper shapes 
of the valence molecular four-component spinors must be restored in atomic core regions after the RECP calculation 
of that molecule performed. 

In 1959, a nonrelativistic procedure of restoration of the orbitals from smoothed Phillips-Kleinman pseudoorbitals 
was proposed l55| based on the orthogonalization of the latters to the original atomic core orbitals. In 1985, Pacios 
& Christiansen [56j suggested a modified orthogonalization scheme in the case of shape-consistent pseudospinors. 
At the same time, a simple procedure of one-center restoration employing the idea of generation of equivalent basis 
sets in four-component Dirac-Fock and two-component RECP/SCF calculations was proposed in [21| (i.e. NOCR 
procedure, see below) and first applied to evaluation of the P,T-odd spin-rotational Hamiltonian parameters in the 
PbF molecule. In 1994, similar procedure was used by Blochl inside the augmentation regions [51J in solids to 
construct the transformation operator between pseudoorbitals ("PS") and original orbitals ("AE") in his projector 
augmented-wave method. 

All the above restoration schemes can be called by "nonvariational" as compared to the "variational" one-center 
restoration (VOCR, see below) procedure proposed in 0, [5^. Proper behavior of the molecular orbitals (four- 
component spinors in the relativistic case) in atomic cores of molecules can be restored in the scope of a variational 
procedure if the molecular pseudoorbitals (two-component pseudospinors) match correctly the original orbitals (large 
components of bispinors) in the valence region after the molecular RECP calculation. This condition is rather correct 
when the shape-consistent RECP is involved to the molecular calculation with explicitly treated outermost core 
orbitals and, especially, when the GRECP operator is used as is demonstrated in 0. 

At the restoration stage, a one-center expansion on the spherical harmonics with numerical radial parts is most 
appropriate both for orbitals (spinors) and for the description of "external" interactions with respect to the core regions 
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of a considered molecule. In the scope of the discussed two-step methods of the electronic structure calculation of 
a molecule, finite nucleus models and quantum electrodynamic terms including, in particular, two-electron Breit 
interaction may be used without problems. 

One-center expansion had been applied first to whole molecules by Desclaux & Pyykko in relativistic and nonrel- 
ativistic Hartree-Fock calculations for the series CH4 to PbEU and then in the Dirac-Fock calculations of CuH, 
AgH and AuH 61] and other molecules [62(. A large bond length contraction due to the relativistic effects has been 
estimated. However, the accuracy of such calculations is limited in practice because the orbitals of the hydrogen atom 
are reexpanded on a heavy nucleus in all the coordinate space. It is worth to note that the RECP and one-center 
expansion approaches were considered earlier as alternatives to each other \(y\ . 

The applicability of the proposed two-step algorithms for calculation of wave functions of molecules with heavy atoms 
is a consequence of the fact that the valence and core electrons may be considered as two subsystems, interaction 
between which is described mainly by some integrated and not by detailed properties of these subsystems. The 
methods for consequent calculation of the valence and core parts of electronic structure of molecules give us a way to 
combine relative simplicity and accessibility of both molecular RECP calculations in gaussian basis set and relativistic 
finite-difference one-center calculations inside a sphere with the atomic core radius. 

The first two-step calculations of the P,T-odd spin-rotational Hamiltonian parameters were performed for the PbF 
radical about 20 years ago [2 lL l65| with the semiempirical accounting for the spin-orbit interaction. Before, only 
nonrclativistic SCF calculation of the TIF molecule using the relativistic scaling was carried out in which 

the P,T-odd values were almost three times underestimated as compared to the relativistic DF calculations. The 
latter were first performed only in 1997 by Laerdahl et al. .68] and by Parpia 0|. The next two-step calculation, for 
PbF and HgF molecules 0], was carried out with the spin-orbit RECP part taken into account using the method 
suggested in ]. 

Later we performed correlation GRECP/NOCR calculations of the core properties in YbF 0], BaF 0], again in 
YbF 74] and in TIF 0]. In 1998, first all-electron Dirac-Fock calculations of the YbF molecule were also performed 
by Quiney et al. 0] and by Parpia [j^. Very recently we finished extensive two-step calculations of the P,T-odd 
properties and hyperfine structure of the excited states of the PbO molecule 0, . 

In the paper, the main features of the used two-step method arc presented and only the last series of the two- 
step calculations are discussed, in which electron correlations are taken into account by a combined method of the 
second-order perturbation theory (PT2) and configuration interaction (CI), or "PT2/CI" [79( (for BaF and YbF), by 
the relativistic coupled cluster (RCC) method jsfl 0] (f° r TIF and PbO), and by the spin-orbit direct-CI method 
0, 13 (f° r PbO). In the discussed ab initio calculations the best up to-date accuracy was attained for the 
hyperfine constants and P,T-odd parameters regarding the molecules containing heavy atoms. 



I. GENERALIZED RELATIVISTIC EFFECTIVE CORE POTENTIAL METHOD. 



When core electrons of a heavy-atom molecule do not play an active role, the effective Hamiltonian with RECP 
can be presented in the form 

H Ef = ]>> Schr M + u Ef (*„)] + £ - 1 - • (i) 

- ^ ■ 'iv iv 

Hamiltonian is written only for a valence subspace of electrons which are treated explicitly and denoted by indices 
i v and j v . In practice, this subspace is often extended by inclusion of some outermost core shells for better accuracy 
but we will consider them as the valence shells below if these outermost core and valence shells are not treated 
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using different approximations. U Ef is an RECP (or relativistic pseudopotential) operator that is usualiy written 
in the separable (e.g., see j^l an d references) or radially- local (semi- local) approximations when the valence 
pseudospinors are smoothed in the heavy-atom cores. The RECP operator simulates, in particular, interactions of the 
explicitly treated electrons with those which are excluded from the RECP calculations. Besides, the generalized RECP 
operator 0,0] described below can be used that includes the radially-local, separable and Huzinaga-type relativistic 
pseudopotentials as special cases. Additionally, the GRECP operator can include terms of other types which are 
important first of all for economical but precise treatment of transition metals, lanthanides and actinides. With 
these terms, accuracy provided by GRECPs can be even higher than the accuracy of the frozen core approximation 
employing the same separation on subspaces of valence and core electrons 

00- In E q- CD, hSchr is the one-electron 

Schrodingcr Hamiltonian 

hSchr = _I^2_^£ (2) 

2 r w 

where Zi C is the charge of the nucleus decreased by the number of inner core electrons. Contrary to the four-component 
wave function used in Dirac-Coulomb(-Breit) calculations, the pseudo-wave function in the (G)RECP case can be both 
two- and one-component. The use of the effective Hamiltonian (|TJ instead of the all-electron relativistic Hamiltonians 
leads to the question about its accuracy. It was shown both theoretically and in many calculations (see and 
references) that the usual accuracy of the radially- local RECP versions is within 1000-3000 cm -1 for transition 
energies between low-lying states. 

The GRECP concept is introduced and developed in a series of papers (see and references) . In contrast 

to other RECP methods, GRECP employs the idea of separating the space around a heavy atom into three regions: 
inner core, outer core and valence, which are treated differently. It allows one to attain practically any desired 
accuracy, while requiring moderate computational efforts. 

The main steps of the scheme of generating the GRECP version with the separable correction taken into account 
are: 

1. The numerical all-electron relativistic calculation of a generator state is carried out for an atom under consideration. 
For this purpose, we use the atomic hfdb code [ssjlsq . 

2. The numerical pseudospinors f n ij(f) are constructed of the large components fnlj(f) of the outer core (OC) and 
valence (V) HFD(B) spinors so that the innermost pseudospinors of them (for each I and j) are nodeless, the next 
pseudospinors have one node, and so forth. These pseudospinors satisfy the following conditions: 

fn»(r) = { fnlj(r) > 5 r ~ R - (3) 



Z = 0,1,...,L, i = |i±||, 

71 71 OCl Tl QC7 . . . , Tb v , 

where n v ,n OCl n' oc are principal quantum numbers of the valence and outer core spinors, L is one more than the 
highest orbital angular momentum of the inner core (IC) spinors. The leading power 7 in the polynomial is typically 
chosen to be close to L+l in order to ensure sufficient ejection of the valence and outer core electrons from the inner 
core region. The coefficients are determined by the following requirements: 

*{fnlj} set is orthonormalized, 

• y and its first four derivatives match / n y and its derivatives, 
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• y is a smooth and nodeless function, and 

• fnij ensures a sufficiently smooth shape of the corresponding potential. 

R c is chosen near the extremum of the spinor so that the corresponding pseudospinor has the defined above number 
of nodes. In practice, the R c radii for the different spinors should be chosen close to each other to generate smooth 
potentials. 

3. The U n ij potentials are derived for each ^=0, . . . , L and j=\l± \ \ for the valence and outer core pseudospinors so that 
the f n ij are solutions of the nonrelativistic-type Hartree-Fock equations in the jj-coupling scheme for a "pseudoatom" 
with the removed inner core electrons. 



UniAr) = f-i](r) 



1 d 2 1(1 + 1) Z ic 
2^2 2r 2 + ~7~ 



J(r) + K(r)+£„ J/ n y(r) - 

E £n'nljfn'lj(r) , (4) 

where J and K are the Coulomb and exchange operators calculated with the f n ij pseudospinors, £„y are the one- 
electron energies of the corresponding spinors, and e n >nij are off-diagonal Lagrange multipliers (which are, in general, 
slightly different for the original bispinors and pseudospinors). 

4. The GRECP operator with the separable correction written in the spinor representation is as 

L 1+1/2 

U Ef = U nvLJ (r) + E E { [UnMr) ~ E>».w(r)]P« 

J=0j=|{-1/2| 

+ E^-yM " U nJj(r)]P noa ,ij +£P nM ,«[t^.,w(r) - U nM {r)] (5) 

n c Hoc 

U noc ,ij(r) + U <c! ij(r) 



n oc ,'n , oc 



where 



Pij = E \ l j m )(lj m \> Pn oc ,lj = E \ n oc,ljm)(n oc ,ljm\, 

m——j m——j 

\ljTn)(ljm\ is the projector on the two-component spin-angular function Xljmi \ n oc^j m ){ n oc^j m \ is the projector 
on the outer core pseudospinors fn oc ,ijXijm-) 

and J = L + 1/2. 

Two of the major features of the GRECP version with the separable correction described here are generating of 
the effective potential components for the pseudospinors which may have nodes, and addition of non-local separable 
terms with projectors on the outer core pseudospinors (the second and third lines in Eq. (J5J) to the standard semi- 
local RECP operator (the first line in Eq. JSJ). Some other GRECP versions are described and discussed in papers 
00,00 in details. 

The GRECP operator in spinor representation J5J is mainly used in our atomic calculations. The spin-orbit 
representation of this operator which can be found in 



is more efficient in practice being applied to molecular 
calculations. Despite the complexity of expression (JSJ for the GRECP operator, the calculation of its one-electron 
integrals is not notably more expensive than that for the case of the conventional radially- local RECP operator. 
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II. NONVARIATIONAL ONE-CENTER RESTORATION METHOD 

The electronic density evaluated from the two-component pseudo-wave function (that, in turn, is obtained in 
the (G)RECP calculation accounting for correlation) very accurately reproduces the corresponding all-electron four- 
component density in the valence and outer core regions. In the inner core region, the pseudospinors are smoothed, 
so that the electronic density with the pseudo-wave function is not correct. When operators describing properties of 
interest are heavily concentrated near the nucleus, their mean values are strongly affected by the wave function in the 
inner region. The four-component molecular spinors must, therefore, be restored in the heavy-atom cores. 

All molecular spinors <p p can be restored as one-center expansions in the cores using the nonvariational one-center 



restoration (NOCR) scheme |2ll|^|7^l7J,l7J,l7J,lll|7J,|7g| that consists of the following steps: 

• Generation of equivalent basis sets of one-ccntcr four-component spinors [ ^ ra ] and of smoothed 

V 9nlj(r)X2j-l,jm J 

two-component pseudospinors fnlj( r )xijm in finite-difference all-electron Dirac-Fock(-Breit) and GRECP/SCF cal- 
culations of the same configurations of a considered atom and its ions. The nucleus is usually modeled as a uni- 
form charge distribution within a sphere. The all-electron four-component hfdb [sfil Isfii Is?! ] and two-component 
GRECP / hfj [g^, I^J codes are employed by us to generate two equivalent numerical basis sets used at the restora- 
tion. These sets, describing mainly the atomic core region, are generated independently of the basis set used for the 
molecular GRECP calculations. 

• The molecular pseudospinorbitals are then expanded in the basis set of one-center two-component atomic pseu- 
dospinors (for r<i?noc n where i? noC r is the radius of restoration that should be sufficiently large for calculating core 
properties with required accuracy), 

Lmax J = l*+1/2| 

M X ) ~ E E J2^nl 3 mfnlAr)Xhm , (6) 
(=0 j = \l-l/2\n,m 

where x denotes spatial and spin variables. Note that for linear molecules only one value of m survives in the sum 
for each (j> p . 

• Finally, the atomic two-component pseudospinors in the molecular basis are replaced by equivalent four-component 
spinors and the expansion coefficients from Eq. ||SJ) are preserved: 

imaxJ = l' + l/2| ( f I \ \ 

k(*)«E E E<U ■ (7) 

1=0 j = \l-l/2\n,m \ 9nlj\ r )X2j-l,jm ) 

The molecular four-component spinors constructed this way are orthogonal to the inner core spinors of the heavy 
atom, because the atomic basis functions used in Eq. Q) are generated with the inner core electrons treated as frozen. 



III. VARIATIONAL ONE-CENTER RESTORATION 



In the variational technique of the one-center restoration (VOCR) [58l |59|, the proper behavior of the four- 
component molecular spinors in the core regions of heavy atoms can be restored as an expansion on the spherical 
harmonics inside the sphere with a restoration radius, R VOC r, that should not be smaller than the matching radius, 
R c , used at the RECP generation. The outer parts of spinors are treated as frozen after the RECP calculation of a 
considered molecule. This method enables one to combine the advantages of two well-developed approaches, molecular 
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RECP calculation on gaussian basis set and atomic-type one-center calculation on numerical basis functions, by the 
most optimal way. It is considered theoretically in |59j and some results concerning the efficiency of the one-center 
reexpansion of orbitals on another atom can be found in [47) . 

The VOCR scheme can be used for constructing the core parts of the molecular spinors (orbitals) instead of using 
equivalent basis sets as is in the NOCR technique. Some disadvantage of the NOCR scheme is that molecular 
pseudoorbitals are usually computed in a spin-averaged GRECP/SCF molecular calculation (i.e. without accounting 
for the effective spin-orbit interaction) and the reexpansion of molecular pseudospinorbitals on one-center atomic 
pseudospinors is yet restricted by accuracy, as it was seen in our GRECP/RCC/NOCR calculations Q of the TIF 
molecule (see below). With the restored molecular bispinors, the two-electron integrals on molecular bispinors can 
be easily calculated when following the scheme presented in Eqs. l|9"))-l|12|) of the next section (see 59] for more 
details). Thus, the four-component transfomation from the atomic basis that is now the most time-consuming stage 
of four-component calculations of heavy-atom molecules can be avoided. 

However, the most interesting direction in development of the two-step method is the possibility to "split" the 
correlation structure calculation of a molecule on two consequent correlation calculations: first, in the valence region, 
when the outer core and valence electrons are explicitly involved into the GRECP calculation, and then, in the core 
region, when the outer and inner core space regions are only treated at the restoration stage. Correlation effects in 
the valence and outer core regions (not only valence parts of molecular orbitals but also configuration coefficients) 
can be calculated, for example, by the GRECP/RCC method with very high accuracy. Then correlation effects in 
the inner and outer core regions can be taken into account using the Dirac-Coulomb(-Breit) Hamiltonian and the 
one-center expansion. Increasing the one-center restoration radius R VOC r one can take into account correlation effects 
in the intermediate region (outer core in our case) with the required accuracy. Roughly speaking, the computational 
efforts for the correlation structure calculations in the core and valence regions are added when using the two-step 
approach, whereas in the conventional one-step scheme, they have multiplicative nature. 



To evaluate one-electron core properties (hyperfine structure, P,T-odd effects etc.) employing the above restoraton 
schemes it is sufficient to obtain the one-particle density matrix, {D pq }, after the molecular RECP calculation on the 
basis of pseudospinors {4> p }. At the same time, the matrix elements {Wp 9 } of the property operator W(x) should 
be calculated on the basis of equivalent four-component spinors {4> P }- The mean value for this operator can be then 
evaluated as 



The only explicitly treated valence shells are used for evaluating a core property when applying Eq. (jSJ) since the atomic 
frozen core (closed) shells do not usually contibute to the properties of practical interest. However, correlations with 
the core electrons explicitly excluded from the RECP calculation can be also taken into account if the effective operator 
technique [gel is applied to calculate the property matrix elements {WE^} on the basis set of bispinors {</> p }. Then, in 
expression (JSJ one should only replace {Wp g } by {WEJ}. Alternatively, the correlations with the inner core electrons 
can be, in principle, considered within the variational restoration scheme for electronic structure in the heavy atom 
cores. Strictly speaking, the use of the effective operators is correct when the molecular calculation is carried out with 
the "correlated" GRECP (see in which the same correlations with the excluded core electrons are taken into 

account at the GRECP generation as they are in constructing {WE?}. Nevertheless, the use of the (G)RECP that 



IV. TWO-STEP CALCULATION OF MOLECULAR PROPERTIES. 




(8) 



p'-i 
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does not account for the core correlations at the molecular calculation stage can be justified if they do not influence 
dramatically on the electronic structure in the valence region. The latter approximation was applied in calculations 
of the YbF and BaF molecules described in the following section. 

When contributions to (W) are important both in the core and valence regions, the scheme for evaluating the mean 
value of W(x) can be as follows: 

• calculation of matrix elements on the molecular pseudospinorbitals (which are usually linear combination of atomic 
gaussians) over all the space region using conventional codes for molecular RECP calculations (though, the operator 
XV can be presented in different forms in calculations with the RECP and Dirac-Coulomb(-Brcit) Hamiltonians), 

W = E 5 ^ / W M to dx ; (9) 



r<oo 



• calculation of the inner part of the matrix element of the operator with the same molecular pseudospinorbitals using 
the one-center expansion for {<fi p } (i? cr stays for i? noC r or i? V ocr below, R OCI > R c ): 

= / to w w to dx ; ( io ) 



™ r<R ocr 

• calculation of the inner part of the matrix element of the operator with the restored molecular four-component 
spinors using the one-center expansion for {<fi p }: 

W < = E 5 ^ / ^<(x) W(x)0<(x)dx. (11) 

Pq r<i? ocr 

The matrix element (W) of the W(x) operator is evaluated as 

(W) = (W) - (Wf + (W) K . (12) 

Obviously, that the one-center basis functions can be numerical (finite-difference) for better flexibility. 

The mean values of the majority of operators for the valence properties can be calculated with good accuracy 
without accounting for the inner parts of the matrix elements (11011 and (flip . As is discussed above, for calculating the 
mean values of the operators singular near nuclei it is sufficient to account only for the inner parts, l|ll(l . of the matrix 
elements of W(x) on the restored functions (j) p (x), whereas the other contributions, © and QlOf) . can be neglected. 



Calculation of properties using the finite-field method |9jJ, |92J and Eq. (JSJ within the coupled-cluster approach is 
described in section IVll 



V. CALCULATION OF SPIN-ROTATIONAL HAMILTONIAN PARAMETERS FOR YBF AND BAF 

MOLECULES 

Effective Hamiltonian. The molecular spin-rotational degrees of freedom of YbF and BaF for the 171 Yb and 137 Ba 
isotopes with nuclear spin i=~ are described by Hamiltonian [|2,lzi3>0] 

H sr = BN 2 + 7SN - D e XE + SAI 

+ W A k A \xS-i+{W s k s + W d d c )S\ , (13) 

where N is the rotational angular momentum, B is the rotational constant, S is the effective spin of the electron 
[9^, I is the spin of the 171 Yb ( 137 Ba) nucleus, A is the unit vector directed along the molecular axis from the heavy 
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nucleus to fluorine. The spin-doubling constant 7 characterizes the spin-rotational interaction. D e and E are the 
molecular dipole moment and the external electric field. The axial tensor A describes magnetic hyperfine structure 
on the ytterbium (barium) nucleus: 

fi N I ■ axr ( a\ 

H h f s = — 3 — , a = , (14) 



which can be also expressed in the notations of the isotropic A=(A\\+2A±) /3 and dipole yLj=(A|| — A±)/3 hyperfine 
structure constants. The hyperfine interaction with the 19 F nucleus is significantly smaller and is not of interest for 
the considered effects. 

The last three terms in Eq. (|13fl account for the P- and P,T-odd effects: the first of them describes the electromag- 
netic interaction of the electron with the anapole moment of the nucleus with the constant kx |3 ! second term 
accounts for the scalar P,T-odd electron- nucleon interaction characterized by the dimensionless constant kg; the third 
term corresponds to interaction of the electron EDM d D with the internal molecular field E mo \: 




H d = 2d c ^ Q ^. E mol , (15) 
W d d e = 2( 2 £ 1/2 | Hd(t)\ 2 Z 1/2 ) . (16) 

i 

The constant \Wd characterizes an effective electric field on the unpaired electron. All the P- and P,T-odd constants 
Wx depend on the electron spin density in the vicinity of the heavy nucleus. The reliability of their calculation can 
be tested by comparison of the calculated and experimental values for the hyperfine constants A and Ad- 
Details of calculations. Let us consider the scheme of calculation for the above molecules on example of the YbF 
calculation carried out mainly in paper |74|. The stage of the GRECP calculation of the valence structure for the 
ground state 2 £ of YbF was performed by analogy with the previous YbF calculation [72J . The main difference from 
the latter was the use of technique |46L l84| for freezing the 5s and 5p pseudospinors derived from the calculation 
on the Yb 2+ cation, ft was necessary to freeze these shells because their polarization was taken into account at the 
stage of calculating the effective operator (EO, see f° r details). The spin-averaged GRECP calculations in 

the framework of the RASSCF method (9^, Ih3 with 5284 configurations were performed for 11 electrons distributed 
in RAS-1=(2,0,0,0), RAS-2=(2, 1,1,0) and RAS-3=(6,4,4,2) active orbital subspaces according to the A ll B ll B 2 and 
A 2 irreducible representations (irreps) of the Ci v symmetry group used in the calculation (see paper 96] for details 
on the RASSCF method). Similar calculation was also carried out for the BaF molecule |73j. 

Expressions for the Wd parameter which characterizes the internal molecular field are presented in Eqs. (|15fl and 
l)16p. whereas expressions for other electronic matrix elements which correspond to the parameters A, A^ and Wx 
of operator l|13|) can be found in papers 

HQ 

All the radial integrals and atomic four-component spinors were 
calculated for the finite nucleus 171 Yb using a model of a uniform charge distribution within a sphere. As is mentioned 
in the introduction, atomic matrix elements of operator (|15f) are proportional to Z 3 . The same scaling is applicable to 
the constant Ws, while the matrix elements which contribute to the constants Wx are proportional to Z 2 . The nuclear 
charge of fluorine is eight times smaller than that of ytterbium, therefore, the contributions to the Wx parameters 
from the vicinity of the fluorine nucleus can be neglected. The additional argument in favour of this approximation 
is connected with the circumstance that the state of the unpaired electron in the YbF molecule is described mainly 
by the ytterbium orbitals and thus the spin density is also localized near ytterbium. 

Results. Results of the calculation for the spin- rotational Hamiltonian parameters are presented in Table HJ Com- 
parison of the results of the GRECP/RASSCF calculations [JJ with the results of the GRECP/RASSCF/EO calcu- 
lations [3 confirms the conclusion made in paper th&t the core- valence correlations of the unpaired YbF electron 
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occupying mainly the hybridized 6s — 6p orbital of Yb with the 5s and then with bp shells of the Yb core (first of all, 
spin polarization of the 5s, bp shells) contribute very essentially to the hyperfine and the P- and P,T-odd constants. 

Our final results for the hyperfine structure constant A differ by less than 3% from the experimental value 0|. This 
means that the space-isotropic part of the unpaired spin density in the vicinity of the Yb nucleus in our calculations 
is obtained with rather good accuracy. However, it is also important to reproduce the anisotropic (dipole) part of the 
spin density which contributes to the dipole constant Ad- The value for Ad obtained in the RASSCF/EO calculation 
is in slightly better agreement with the experimental datum than the value from paper [7^ but is still underestimated 
by 23%. About one half of this difference can be explained by the fact that the 4/ shell of Yb was frozen in our 
molecular calculation (it was first pointed out by Khriplovich that accounting for excitations from the 4/-she ll can 
be important). In particular, these excitations can explain the small value of the spin-doubling constant 7 |9^.ll0c| . 
It was shown in semiempirical calculation [9^] that the contribution of the 4/-shell excitations to the spin density 
can result in a significant correction to the constant Ad- Using Eqs. (19) and (31) from paper j^5], one can go to the 
following estimates for the contributions from the 4/-shell excitations to A and A&: 

SA « -3 MHz, 5A d rs 15 MHz. (17) 

Note that the semiempirical corrections obtained in [j^ arise from the admixture to the molecular wavefunction of 
the configuration with the hole in the 4/-shell. The weight of this configuration was estimated in paper [9^ to be of 
the order of 4%. This admixture is a purely molecular effect and is not described by the effective operator technique. 
Thus, one can conclude that the 5 Ad value can be added to the RASSCF/EO value for Ad, resulting in Ad ~ 94 MHz 
which is essentially closer to the experimental value of 102 MHz. Note also that ab initio correlation calculation 
of YbF with, at least, 4/ shell treated explicitly should be performed to check the above described mechanism of 
increasing the calculated Ad value. It is likely that the real situation is more complicated and indirect contributions 
to Ad caused by the 4/ shell correlation-relaxation effects and by other higher order contributions (which are not 
accounted for in the considered calculation) are mainly responsible for the obtained underestimation of Ad- 
it is essential that similar contribution of the 4/-shell excitations to the constant Wd is strongly suppressed. Indeed, 
operator l|15l) mixes /- and d-waves. The 4d-shell is localized rather deeply in the ytterbium core and its mixture 
with the 4/-orbitals by the molecular field is very small while the bd electrons are relatively weakly-bound and does 
not penetrate essentially into the region of the 4/-shell localization. Similar contributions to other constants Wx are 
negligible due to the contact character of the corresponding interactions. 

One can see that the values of the Wd constant from the unrestricted DHF calculation (accounting for the spin- 
polarization of core shells) (t^i the latest semiempirical calculation 95] and the GRECP/RASSCF/NOCR/EO calcu- 
lation are in very good agreement. It is also important for reliability of the carried out calculations that the contribu- 
tion from the valence electron to Wd in paper [76J differs only by 7.4% from the corresponding GRECP/SCF/NOCR 
calculation jj^ (see Table QJ. Another DHF calculation |75[ gives a value that is about two times smaller (that is 
rather explained by different definitions used for Wd)- 

Similar increase in the values for the hyperfine constants and parameters of the P,T-odd interactions when the 
correlations with the core shells (first of all, 5s, bp) are taken into account is also observed for the BaF molecule 
[73I as one can see in Table . Of course, the corrections on the 4/-electron excitations are not required for this 
molecule. The enhancement factor for the P,T-odd effects in the BaF molecule is three times smaller than in YbF 
mainly because of smaller nuclear charge of Ba. Therefore, barium fluoride can be considered as a less promising 
molecule for a search for the PNC effects. 
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TABLE I: Parameters of the spin-rotational Hamiltonian for the 171 YbF and 137 BaF molecules; A = (Ay +2A±)/3 (isotropic) 
and Ad = (j4.ii — A±)/3 (dipole) are the HFS constants. 





.4 


A d 


w d 


Wa 


W s 




(MHz) 


(MHz) (10 25 gS^ 


) (kHz) 


(Hz) 


The 171 YbF molecule 


Experiment [98J 


7617 


102 








Semiempirical [101] 






-1.5 


730 


-48 


Semiempirical [731 (with 4/-correction) 






-1.26 




-43 


GRECP/SCF/NOCR [72] 


4932 


59 


-0.91 


484 


-33 


GRECP/RASSCF/NOCR [72J 


4854 


60 


-0.91 


486 


-33 


Restricted DHF (Quiney, 1998) [75] 


5918 


35 


-0.31 


163 


-11 


Restricted DHF + core polarization 


7865 


60 


-0.60 


310 


-21 


Rescaled restricted DHF 






-0.62 


326 


-22 


Unrestricted DHF (Parpia, 1998) [76] 












(unpaired valence electron) 






-0.962 






Unrestricted DHF [22] 






-1.203 




-22 


GRECP/RASSCF/NOCR/EO [74] 


7842 


79 


-1.206 


634 




GRECP/RASSCF/NOCR/EO [74] 












(with 4/-correction) 


7839 


94 


-1.206 


634 




The 137 BaF molecule 


Experiment [1021° 


2326 


25 








Semiempirical [20J 






-0.41 


240 


-13 


Experiment [1031'' 


2418 


17 








Semiempirical [20J 






-0.35 


210 


-11 


GRECP/SCF/NOCR [73] 


1457 


11 


-0.230 


111 


-6.1 


GRECP/RASSCF/NOCR [£3J 


1466 


11 


-0.224 


107 


-5.9 


GRECP/SCF/NOCR/EO [73] 


2212 


26 


-0.375 


181 




GRECP/RASSCF/NOCR/EO [73] 


2224 


24 


-0.364 


175 





a The HFS constants are measured in an inert gas matrix 1^21 and the semiempirical Wx values are based on these data. 
'The HFS constants were measured in a molecular beam 103]. 



VI. CALCULATION OF 205 TLF MOLECULE. 



Effective Hamiltonian. Here we consider the interaction of the EDM of unpaired proton in 205 Tl with the internal 

electromagnetic field of the 205 T1F molecule Q|. This molecule is one of the best objects for the proton EDM 
measurements. The effective interaction with the proton EDM in TIF is written in the form 

H cS = (d v + d M )T/I • A , (18) 

where / is the Tl nucle ar spin operator, A is the unit vector along z (from Tl to F), d v and d are volume and 
magnetic constants |l04| 

d v = 6SX = (-d p R + Q)X , (19) 
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S is the nuclear Schiff moment, d p is the proton EDM, 

} 1 

(20) 



9 (A 



- x,y,z—0 

pip(r) is the electronic density calculated from the electronic wave function ip; 

d M = 2V2(d p + d N )(^ + -^-]M , (21) 
\Z 2mc J 

where djv is the nuclear EDM arising due to P,T-odd nuclear forces; fj,, m and Z are the magnetic moment, mass and 
charge of the Tl nucleus; c is the velocity of light, 

M=^lZ(^) \4>), (22) 

li is the orbital momentum for electron i; cti are its Dirac matrices. R and Q are factors determined by the nuclear 
structure of 205 T1: 

R = (V>iv(f„)| YM«I Z - S n,3s)r 2 n \^ N (fn)) , (23) 
n 

Q = [S/5(i) N (f n )\} (q n f n )\-ij} N (f n )) - 



l/Z{i/, N (f»)\ Yi<lnrl)\i>N{?n)) (24) 



n 

x ( ipN (r„ ) | V (q n f n ) /rl \ (f„ ) ) ] f 



where ^jv(r n ) is the nuclear wave function. 

Accounting for H e g leads to different hyperfine splitting of TIF in parallel and antiparallel electric and magnetic 
fields. The level shift hv = A(d v + d M )(I/I ■ A) is measured experimentally (for the latest results on TIF see [l3j)- 

The parameters X of Eq. (|2(J|I and M of Eq. H22fl are determined by the electronic structure of the molecule. They 
were calculated recently for the ground + state of TIF by Parpia j^Sl arid by Laerdahl, Saue, Faegri Jr., and Q uiney 
|6^ . using the Dirac-Fock method with gaussian basis sets of large sizes (see Table ITTf) . (Below we refer to paper Il05 
the LlcuJatio™ „ in de«a ils and not to the brief eonnnnnic.tion 6 8] oT.hc sam e aut hors.) In paper |b| 
the first calculation of 205 TIF that accounts for correlation effects was performed (see also |lod | where the limit on 
the Schiff moment of 205 Tl was recalculated). 



Details of calculations. In paper 14] 21-electron GRECP 107 for Tl was used. The correlation spin-orbital basis 
set consisted of 26s, 25p, 18d, 12/, and lOg gaussian-type orbitals on Tl, contracted to 6s6p4d2flg. The basis was 
optimized in a series of atom ic two-component GRECP calculations, with correlation included by the relativistic 
coupled cluster method |l0s| with single and double cluster amplitudes; the average e nergy of t he two lowest states 
of the atom was minimized. The basis set generation procedure is described in Refs. [lod llicj . The basis set was 
designed to describe correlation in the outer core 5s and 5p shells of Tl, in addition to the 5d and valence shells. 
While correlating the 5s and 5p shells may not be important for a majority of chemical and physical properties of 
the atom, it is essential for describing properties coming from the inner Tl region, including P,T-odd effects. The 
(14s9p4cZ3/)/[4s3p2dl/] basis set from the ANO-L library [97J is used for fluorine. 

First, one-component SCF calculation of the (la . . . 7cr) 14 (l7r27r37r) 12 (l(5) 4 ground state of TIF is performed, using 
the spin-averaged GRECP for Tl which simulates the interactions of the valence and outer core (5s5p5d) electrons 
with the inner core [Kr]4d^ 2 4d^ 2 4/|^ 2 4/^ 2 . This is followed by two-component RCC calculations, with only single 
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(RCC-S) or with single and double (RCC-SD) cluster amplitudes. The RCC-S calculations with the spin-dependent 
GRECP operator take into account effects of spin-orbit interaction at the level of the one-configurational SCF-type 
method. The RCC-SD calculations include, in addition, the most important electron correlation effects. 

At the restoration stage, the nucleus was modeled as a unifor m cha rge distribution within a sphere with radius 
7"nuci = 7.1 fm = 1.34 x 10~ 4 a.u., whereas previous calculations [oil llflfJ ] employed a spherical gaussian nuclear charge 
distributi on (t he root mean square radius in all calculations is 5.5 fm, in accord with the para metrization of Johnson 
and Soff |lll| . and agrees with the experimental value 5.483 fm for the 205 T1 nucleus |ll2|). The two equivalent 
[15sl2pl2e?8/] numerical basis sets were generated for the restoration. 

The quality of the approximation for the two-center molecular spinors and, consequently, of the calculated properties 
increases with the value of £ max . A series of calculations of the M parameter was performed using Eq. Q with basis 
functions including up to p, d and / harmonics. It is found (see Table ITT|) that accounting for only s and p functions 
in the expansion determines M with 90% accuracy. Because the contribution of / functions is only about 0.3% and 
amplitudes of higher harmonics on the nucleus are suppressed by the leading term ~ 7-0' _1 / 2 ) ) the error due to the 
neglect of spherical harmo nics beyond / is estimated to be below 0.1%. Calculation of the X parameter requires s 
and p harmonics (see Ref. |l05^ . although, strictly speaking, d harmonics also give nonzero contributions. 

In the RCC calculations, the X and M parameters were calculated by the finite field method (e.g., see Refs. joiUo^ ). 
The operator corresponding to a desired property (Eqs. (|20[) and iJUJl) is multiplied by a small parameter A and added 
to the Hamiltonian. The derivative of the energy with respect to A gives the computed property. This is strictly 
correct only at the limit of vanishing A, but it is usually possible to find a range of A values where the energy is almost 
linear with A and energy changes are large enough to allow sufficient precision. The quadratic dependence on A is 
eliminated in the present calculations by averaging absolute energy changes obtained with A of opposite signs. 

Results. Calculations were carried out at two internuclear separations, the equilibrium R e — 2.0844 A as in 
Ref. |sa, and 2.1 A, for comparison with Ref. |l05| . The results are collected in Table [n] The first point to 
notice is the difference between spin-averaged SCF values and RCC-S values, the latter include spin-orbit interaction 
effects. These effects increase X by 9% and decrease M by 21%. The RCC-S function can be written as a single 
determinant, and results may therefore be compared with DHF values, even though the RCC-S function is not 
variational. GRECP /RCC-S/N OCR values of the M parameter are indeed within 3% and 1% of difference with the 
corresponding DHF values p^L Il05j (Table ITI]) . This agreement confirms the validity of the used approximations. 
In particular, freezing the inner core shells is justified, as inner core r elaxa tion effects have little influence on the 
properties calculated here, a conclusion already drawn by Quiney et al. |l05| . 

Much larger differences occur for the X parameter. There are also large differences between the two DHF calcu- 
lations for X, which cannot be explained by the small change in internuclear separation. The value of X may be 
expected to be less stable than M, bec ause it is determined by the derivative of the electronic density at the Tl 
nucleus and involves large cancellations 105] between contributions of large and small components, each of them is 
about 20 times larger than their sum. Thus, a strong dependence of X on the basis used may be ex pected. The DHF 
values collected in Table ITT1 indeed show such dependence. Results obtained in Refs. 69] and |l05j with comparable 
even-tempered basis sets, (28s2 8pl2d8/) and (28s28pl4d8/), are rather close, differing by 340 a.u. Improving the 
Tl basis to (34s34pl6<i9/) 105] increases X by 650 a.u. or 8%. Further improvement of the basis may be expected 
to yield even higher X values. The numerical basis functions obtained in atomic DHF calculations and used for the 
restoration are highly acc urat e near the nucleus, so that the GRECP/RCC-S/NOCR value for X, which is h igher 
than that of Quiney et al. |l05| . seems reasonable. The different nuclear models used in the present and DHF |6gL Il05j 
calculations may also contribute to the disagreement in X, which is determined by the derivative of the electronic 
charge density at a single point, the Tl origin. M is affected by t/j in a broader region, and is therefore far less sensitive 
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TABLE II: Calculated X and M JUJl parameters (in a.u.) for the 205 T1F ground state, compared with DHF values |6allOI 
GRECP/RCC-S results include spin-orbit interaction, and GRECP/RCC-SD values also account for electron correlation. 





R e = 2.0844 A 


R = 2.1 A 


Expansion 


s,p 


s,p,d s,p,d,f 


s,p 


s,p 


s,p,d,f 


s, p 








M 


X 




M 


X 


SCF (spin-averaged) 


19.67 


17.56 17.51 


8967 


19.52 


17.43 


8897 


GRECP/RCC-S 


16.12 


13.84 


9813 


16.02 


13.82 


9726 


DHF [69] 


Tl:(28s28pl2d8/) 


15.61 


7743 






DHF [105] 


Tl:(25s25pl2d8/) 










13.64° 


8098 




Tl:(28s28pl4d8/) 










13.62° 


8089 




Tl:(31s31pl5d8/) 










13.66° 


8492 




Tl:(34s34pl6d9/) 










13.63° 


8747 


GRECP/RCC-SD 


11.50 


7635 







a M is calculated in Ref. 105] using two-center molecular spinors, corresponding to infinite L max 111 Eq. Q. 



to the nuclear model. 

The electron correlation effects are calculated by the RCC-SD method at the molecular equilibrium internuclear 
distance R e . A major correlation contribution is observed, decreasing M by 17% and X by 22%. 

The hyperfine structure constants of Tl &p\/ 2 an d Tl 2+ 6s 1 , which (like X and M) depend on operators concentrated 
near the Tl nucleus, were also calculated. The errors in the DF values are 10-15%; the RCC-SD results are within 
1-4% of experiment. The improvement in X and M upon inclusion of correlation is expected to be similar. 



VII. CALCULATIONS OF 2(J7 PBO MOLECULE. 

It was noted previously that the experiments on the excited a(l) Q or B(l) |l 13| states of PbO having nonzero 
projection of total electronic momentum on the internuclear axis can be, in principle, sensitive enough to detect d e 
three or even four orders of magnitude lower than the current limit. The knowledge of the effective electric field, Wd, 
seen by an unpaired electron is required for extracting d e from the measurements. In p aper s | t7L It8I | . W d for the a(l) 
and B(l) states of the PbO molecule was calculated by the RCC-SD 0E3 and CI 114, 82, 3/3 methods. To 
provide an accuracy check in calculation of the electronic structure near the Pb nucleus the hyperfine constant, Am , 
was also calculated. 



Details of calculations. The 22-electron GRECP for Pb |107j is used at the first stage of the two-step calculations 
of PbO: the inner shells of the Pb atom (Is to 4/) are absorbed into the GRECP, and the 5s5p5d6s6p electrons 
and all the oxygen electrons are treated explicitly at the integral preparation part. Two series of RCC-SD calcu- 
lations with 10 (when freezing the 5s5p5d shells of lead employing the level shift technique 0] and freezing the 
Is shell of oxigen after the RASSCF stage) and all 30 electrons treated explicitly were performed. Only 10 elec- 
trons were treated explicitly in the CI calculations. For each series of calculations, correlation spin-orbital basis 
sets are optimized in atomic two-component GRECP/RCC calculations of Pb. The generated basis sets on Pb are 
(14sl8pl6d8/)/[4s7p5d3/] for 30-electron and (15sl6pl2d9/)/[5s7p4d2/] for 10-electron calculations. The correlation- 
consistent (10s5p2dl/)/[4s3p2dl/] basis listed in the MOLCAS 4.1 library 97] is used for oxygen. We found that 
the / basis function of oxygen has little effect on the core properties calculated here and it was not included into the 
RCC calculations to reduce expenses. 

The leading AS coupling terms and configurations for the a(l) and B(l) states are 3 S + a\u\u\'r:\T:\ and 3 IIi 
(T i (T 2 cr 3 7r i 7r 2; correspondingly. The molecular orbitals used in the CI and RCC-SD calculations are obtained by the 
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RASSCF method 0,0 with the spin-averaged GRECP part j^li i- e - on ly scalar-relativistic interactions are taken 
into account in the RASSCF calculation. The a(l) state is of practical interest first of all and, therefore, the lowest 
3 E+ state was calculated at this stage. Using the C 2)J point group classification scheme, two Ai orbitals are in RAS1, 
six orbitals (two Ai, two Bi, and two B 2 ) in RAS2, and 50 (20 Ai, 6 A 2 , 12 Bi, and 12 B 2 ) in RAS3 for 10-electron 
calculation. 2 Ai orbitals are in RAS1, 6 orbitals (2 Ai, 2 Bi, and 2 B 2 ) in RAS2, and 41 (16 Ai, 5 A 2 , 10 B x , and 
10 B 2 ) in RAS3 for 30-electron calculation. No more than two holes in RAS1 and two particles in RAS3 are allowed. 

At the restoration stage the nucleus is modeled by a uniform charge distribution within a sphere of radius r nuc i = 
7.12 fm = 1.35 x 10~ 4 a.u . The root mean square radius of the nucleus is 5.52 fin, in accord with the parametrization 
of Johnson and Soff 111|, and agrees with the experimental value of 5.497 fm 112 for the 207 Pb nucleus. Taking this 
value for the root mean square radius and a Fermi distribution for the nuclear charge changes An and W d by 0.1% or 
less. The equivalent basis sets generated are [9slAp7d] for 10 and [6s7p5d] for 30-electron calculation correspondingly. 

Then the RCC-SD method and the spin-orbit CI approach with the selected single- and double-excitations from 
some multiconfigurational reference states ( "mains" ) |84l Ill4 1 are employed in the sets of different AS many-electron 
spin- and space-symmetry adapted basis functions 

Results. CI calculations |78| were performed at two internuclear distances, R = 3.8 a.u., as is in RCC calculations, 
and R — 4.0 a.u. (in the RCC calculations [73 the only internuclear distance R = 3.8 a.u. was used because of problem 
with convergence at larger distances). The calculated values with the one-center expansion of the molecular spinors 
in the Pb core on either s, s]p or s;p; d partial waves are collected in Table ITTT1 

Let us consider the 5s, 5p, 5d orbitals of lead and Is orbital of oxigen as the outercore ones and a±, cr 2 , 173, tti, 
7r 2 orbitals of the PbO (consisting mainly of 6s, 6p orbitals of Pb and 2s, 2p orbitals of O) as valence. Although in 
the CI calculations we take into account only the correlation between valence electrons, the accuracy attained in the 
CI calculation of A\\ is much better than in the RCC-SD calculation. The main problem with the RCC calculation 
was that the used Fock-space RCC-SD version was not optimal in accounting for the nondynamic correlations (see 
|l09| for details of RCC-SD and CI calculations of the Pb atom), though the potential of the RCC approach for 
electronic structure calculations is very high in prospect, especially, in the framework of the intermediate Hamiltonian 
formulation [si| . 

Then we estimate the contribution from correlations of valence electrons with outercore ones (which also account 
for correlations between outercore pairs of el ectro ns) as difference in the results of the corresponding 10- and 30- 
electron GRECP/RCC calculations (see also |l09j where this correction is applied to the Pb atom). We designate 
such correlations in Table ITTT1 as "outercore correlations". When taking into account outercore contributions at the 
point R = 4.0 a.u. we used the results of the RCC calculation at the point R = 3.8 a.u. Since these contributions 
are relatively small and because the correlations with the outercore electrons are more stable than correlations in the 
valence region when the internuclear distance is changed the used approximation seems us reasonable. At last, the 
linea r extrapolation of the result s of c alculations to the experimental equilibrium distances, R e — 4.06 a.u. for a(l) 
[n3 and R e = 3.91 a.u. for 5(1) [l3 is applied. The final results are: An = -3826 MHz, W d = -6.1-10 24 Hz/(e • cm) 



for a(l) and An = 4887 MHz, W d = -8.0-10 24 Hz/(e • cm) for 5(1) states. The estimated error for the final W d 



Ml "■ " ' ' 

parameter is 20% for B(l) state. For a(l ) the estimated error bounds put the actual W d value between 90% and 

130% of the final value (for details see |78j1. 
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TABLE III: Calculated parameters A\\ (in MHz) and W d (in 10 24 Hz/(e ■ cm)) for the a(l) and B(l) states of 207 PbO in the 
internuclear distances 3.8 and 4.0 a.u. T he ex perimental value of Am in a(l) is —4113 MHz |ll7| ]. The preliminary experimental 
value of Ay for 5(1) is 5000 ± 200 MHz [usf . 



State 

Parameters 


/ n \ 22231 3v 
a(l) <7 1 CT 2 a 3 Tr 1 TC 2 ^1 

A || W d 




B(l) a\<rl<j\ 






Expansion s 


s,p SjPj - s iP s >P>d 


s 


s,p s,p,d 


s,p 


s,p,d 




Internuclear distance R. — 3.8 a.u. 










lOe-RASSCF -894 
30e-RASSCF -759 


-1505 -1503 0.73 0.70 
-1705 -1699 0.96 0.91 




1900 


0.0 
0.0 


0.0 
0.0 


lOc-RCC-SD [77] 
30e-RCC-SD [77] 

outercore (30e-RCC-SD - lOe-RCC-SD) 


-2635 -2.93 -3.05 
-2698 -4.10 
-63 -1.05 




3878 
4081 

203 


-11.10 
-9.10 


-10.10 
-9.70 
0.40 


lOc-CI [£8J 
FINAL [781 

(lOc-CI + outercore) 


-3446 -4.13 

-3509 -5.18 
Internuclear distance R — 4.0 a.u. 




4582 
4785 




-10.64 
-10.24 


lOe-RASSCF -770 


-1384 -1383 1.05 1.00 






0.0 


0.0 


lOc-CI [zsj 

FINAL [78] 

(lOe-CI + outercore) 12 


-3689 -4.81 
-3752 -5.86 




4762 
4965 




-7.18 
-6.78 


"It is assumed that the outercore contribution at the internuclear distance i?=4.0 a.u. is 


approximately the same as 


is at R=3.8 


a.u. 



VIII. CONCLUSIONS 

The P,T-parity nonconservation parameters and hyperfine constants are calculated for the heavy-atom molecules 
which are of primary interest for modern experiments on search for the PNC effects. It is found that accounting for the 
electron correlations is necessary for precise calculation of these properties. The developed two-step (GRECP/NOCR) 
scheme of calculation of the properties heavily concentrated in atomic cores is a very efficient way to take account of 
these correlations with moderate efforts. The results of the two-step calculations for the hyperfine constants differ less 
than 10% from the corresponding experimental data. The comparable level of accuracy is expected for the P,T-odd 
spin-rotational Hamiltonian parameters which can not be obtained experimentally. The precision of the discussed 
calculations is limited by the current possibilities of the correlation methods and codes and not by the GRECP 
and NOCR approximations, despite the GRECP/NOCR method allows one to reduce seriously the expenses of the 
correlation treatment as compared to conventional Dirac-Coulomb(-Breit) approaches when using basis of molecular 
spin-orbitals instead of spinors etc. In turn, the attained accuracy is sufficient for a reliable interpretation of the 
results measured in the molecular PNC experiments. 
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